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ABSTRACT 


The  open-water  performance  and  the  pressure  distribu¬ 
tion  on  the  blades  of  ITTC  Model  Propeller  MP  282  (diameter 
of  950  mm)  operating  in  uniform  flow  were  computed  using  the 
Computer  Code  PUF-2,  originally  developed  at  Massachusetts 
Institute  of  Technology  (MIT).  The  computed  results  are  com¬ 
pared  with  experimental  measurements  made  at  Ishikawa jima- 
Harima  Heavy  Industries  (IHI)  Ship  Model  Basin  in  Japan.  The 
predicted  thrust  and  torque  coefficients  were  about  10  percent 
less  than  the  experimental  values  over  the  range  of  advance 
coefficients.  In  general,  the  predicted  pressure  distributions 
are  in  satisfactory  agreement  with  the  experimental  measurements. 


ADMINISTRATIVE  INFORMATION 

This  work  was  funded  by  the  Propulsor  Technology  Subprogram  of  the  6.2  Ships, 
Subs  and  Boats  Program  and  performed  at  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  (DTNSRDC)  under  Work  Unit  1540-002. 


INTRODUCTION 

In  this  report  the  open-water  performance  and  the  pressure  distribution  on 
the  blades  of  ITTC  Model  Propeller  MP  No.  282  operating  in  uniform  flow  were  com¬ 
puted  using  the  computer  code  PUF-2,  originally  developed  at  MIT.1  The  computa¬ 
tional  algorithm  for  the  pressure  distribution  was  developed  and  implemented  to  the 

2 

existing  code,  PUF-2,  by  Kobayashl.  The  computed  results  were  compared  with  the 

3 

experimental  measurements  made  at  IHI  Ship  Model  Basin. 

COORDINATE  SYSTEM  AND  BLADE  GEOMETRY 

The  blade  outline  and  the  principal  characteristics  of  Model  Propeller  MP  282 
are  shown  in  Figure  1  and  Table  1,  respectively.  The  coordinate  system  used  in  PUF-2 
is  shown  in  Figure  2.  In  Figure  3,  the  definition  of  the  blade  surface  ordinates 
provided  by  the  ITTC  Propeller  Committee  (ITTC  PC)  and  that  of  the  PUF-2  are  compared 
According  to  the  ITTC  PC  definition,  the  chordlength  was  defined  along  the  "profile 
base  line”  and  the  pitch  angle  was  defined  by  the  angle  between  the  profile  base  line 
and  the  line  in  the  plane  of  the  propeller  rotation.  In  PUF-2  coordinate  system,  the 
chordlength  was  defined  along  the  nose-tail  line  and  the  pitch  angle  was  defined  by 
the  angle  between  the  nose-tail  line  and  the  line  in  the  plane  of  propeller  rotation. 


*A  complete  listing  of  references  is  given  on  page  7. 


Therefore,  the  following  coordinate  transformation  was  made  to  convert  ITTC  PC 
data  in  oxy  coordinate  system  into  the  o'x'y'  system  used  in  PUF-2  (see  Figure  3): 


x'  ■  x  cos  a 

-  y 

sin  a 

+  yo 

sin  a 

(1) 

y'  -  x  sin  ct 

+  y 

cos  a 

•  yo 

cos  a 

(2) 

where  yQ  is  the  ordinate  at  the  leading  edge  in  oxy  system,  and  a  is  the  Induced 
pitch  angle  due  to  the  above  coordinate  transformation  defined  by 

y 

a  -  tan  *  —  (3) 

c 

This  coordinate  transformation  results  in  the  new  chordlength,  pitch  angle, 
skew,  and  rake  in  accordance  with  the  PUF-2  definition  (see  Figures  1  and  2): 


■  c/cos  a 

(A) 

*p  +“ 

(5) 

i  PQ  cos  <J>p  -  P'Q  cos(4>p+a) 

(6) 

PQ  sin  4>p  -  P'Q  sin  (4>p-K») 

(7) 

A  small  computer  program  was  developed  to  convert  the  original  data  into  the 
data  for  PUF-2  input*  Table  2  shows  the  radial  distribution  of  the  pitch,  rake, 
skew,  and  the  chordlength  of  the  new  geometry  after  the  coordinate  transformation. 

In  Table  3,  the  original  table  of  offsets  (in  mm)  as  provided  by  ITTC  PC  were  shown. 
Table  4  shows  the  new  table  of  offsets  about  the  nose-tail  line  after  the  coordinate 
transformation.  In  Table  5,  the  same  data  as  in  Table  4  are  shown  together  with 
camber  and  thickness  distribution  at  different  fractions  of  chord  from  Table  4. 

The  17  new  fractions  of  chord,  x£  (0.0,  0.01,  0.025,  0.05,  0.1,  0.2,  0.3,  0.4,  0.5, 
0.6,  0.7,  0.8,  0.9,  0.95,  0.975,  0.99,  1.0)  were  required  by  PUF-2  input. 
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RESULTS  AND  DISCUSSION 

Since  the  original  PUF-2  code  was  written  for  a  propeller  geometry  with  fixed 
thickness  and  camber  distribution  throughout  the  radius  (only  the  maximum  thickness 
and  camber  change  radially),  the  algorithm  had  to  be  modified  slightly  in  order  to 
handle  the  blade  geometry  with  radially  varying  camber  and  thickness  distribution 
such  as  the  present  propeller  (see  Table  5).  The  open-water  performance  was  calcu¬ 
lated  in  terms  of  K^,  K^,  and  efficiency  at  6  different  advance  coefficients: 

J  “  0.9,  0.95,  1.0,  1.054,  1.1,  and  1.163.  In  the  computations  a  value  of  0.007 
was  used  for  the  drag  coefficient  for  all  J  values.  The  computed  values  of  and 
Kq  were  compared  with  experimental  results  in  Table  6.  The  predicted  KT  and 
were  about  15  percent  and  10  percent,  respectively,  lower  than  the  experimental 
values  over  the  range  of  advance  coefficients.  The  efficiency  was  about  7  percent 
lower  than  the  experimental  results.  The  discrepancies  in  open-water  performance 
are  somewhat  larger  than  some  previous  correlation  results.^ 

2 

Some  slight  modifications  were  made  to  the  algorithm  of  Kobayashi  in  order 
to  improve  the  pressure  prediction.  With  this  corrected  version  of  the  program, 
the  pressure  distribution  on  the  same  blade  (Propeller  4498)  as  in  Reference  2 
was  calculated  and  compared  with  other  theory.  The  correlation  was  improved. 

The  pressure  distribution  on  the  blade  of  Propeller  MP  282  was  calculated 
for  a  range  of  J  values.  The  computed  results  were  shown  in  Tables  7  through  12 
for  each  J  value.  The  pressure  coefficients  were  calculated  on  both  the  suction  and 
the  pressure  sides  at  selected  radii  (r/R  -  0.3,  0.4,  0.5,  0.6,  0.7,  0.8,  0.9). 

The  pressure  coefficients  were  nondlmenslonalized  by  the  local  inflow  at  each 
radius : 


where  VR 


Yva2  +  (2*nr)2  . 


(8) 


In  these  tables,  the  lift  coefficients,  C^>  are  also  shown  at  each  radius.  The  lift 
coefficient  was  computed  by  integrating  the  difference  in  pressure  coefficients: 


C 


L 


AC  dx 
p  c 


(9) 
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For  each  radius,  was  computed  at  discrete  fractions  of  chord,  from  2.5  percent 

to  92.5  percent  of  the  chord.  Therefore,  at  the  leading  and  the  trailing  edges, 

AC  was  arbitrarily  set  to  zero  and  the  discrete  values  of  AC  were  curve-fitted 
P  p 

using  cubic-spline  functions  before  integration. 

Table  13  shows  the  experimental  measurements  of  the  pressure  distribution 

3 

made  at  IHI  Ship  Model  Basin.  In  Figures  4  and  5,  the  experimental  measurements 

and  the  predicted  C^  were  compared  at  J  ■  1.054  and  J  »  1.163,  respectively.  The 

experimental  measurements  were  made  at  Reynolds  number,  Rr  *  1.9  x  10^.  The 

calculated  pressure  coefficients  are  in  agreement  with  measurements  on  the  pressure 

side  except  near  the  leading  edge,  but  generally  overpredict  the  suction  side 

pressure.  The  agreement  at  reduced  J  value  is  better  than  that  at  Increased  J  value. 

In  general,  the  predicted  values  are  in  agreement  with  the  experimental  measurements 

throughout  the  radius  at  two  different  J  values. 

In  Figure  6,  the  oil-film  test  results  reconstructed  from  the  photographs  in 

6  6 

Reference  3  are  shown  at  two  different  Reynolds  numbers;  1.1  x  10°  and  2.6  x  10  , 

respectively.  It  is  observed  that  the  suction  side  flow  at  R  ■  1.1  x  10°  is 

6  n 

somewhat  different  from  the  flow  at  R  ■  2.6  x  10  .  At  the  reduced  R  condition, 

n  n 

the  flow  patterns  on  the  suction  side  have  significantly  reduced  shear  stress 
over  the  forward  part  of  the  blade  and  a  clear  separation  occurs  slightly  past 
midchord.  On  the  pressure  side,  reduced  shear  regions  occur  toward  the  leading 
edge  and  some  indication  of  a  leading-edge  laminar  separation  bubble  occurs  at 
both  Reynolds  numbers. 

No  surface  flow  patterns  are  presented  in  Reference  3  for  the  test  R  of 

6  ^ 

1.9  x  10  .  However,  it  is  possible  that  separation  occurred  near  0.7  fraction 

of  chord  on  the  suction  side  and  at  the  leading  edge  on  the  pressure  side  in  the 

form  of  a  bubble.  Such  separation  would  explain  the  suction  peak  on  the  pressure 

side  near  the  leading  edge  and  the  pressure  peak  measured  at  0.7  radius  at  0.7 

fraction  of  chord  (measurements  were  not  made  at  a  similar  chordwlse  position  at 

other  radii)  (see  Figures  4  and  5).  It  is  further  hypothesized  that  the  suction 

side  separation  is  a  thin  layer  with  only  minor  influence  on  the  pressures  away 

from  the  separation  line. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  steady  pressure  distribution  on  a  large  diameter  model  propeller  predicted 
by  a  modification  to  the  Computer  Code  PUF-2  is  in  satisfactory  agreement  with 
experimental  measurements*  Viscous  effects  such  as  suction-side  Separation  and 
leading-edge  laminar  bubble  separation  on  the  pressure  side  may  have  occurred 
during  the  experimental  pressure  measurements.  However,  it  appears  that  the 
viscous  phenomena  have  only  a  local  effect  on  the  pressure  distribution. 

In  order  to  assess  the  validity  of  the  modified  computer  code  for  blade 

pressure  distribution,  it  is  recommended  that  more  comparisons  be  performed  for  a 

4 

wide  range  of  propellers  and  operating  conditions.  Recently,  Greeley  and  Kerwin 
developed  a  more  advanced  numerical  lifting  surface  theory  and  a  computer  code  PSF 
for  the  prediction  of  propeller  steady  performance.  It  is  recommended  that  this 
code  be  modified  to  provide  for  prediction  of  pressure  distribution.  Since  steady 
performance  predictions  by  PSF  were  shown  to  be  in  better  agreement  with  measure¬ 
ments  than  those  by  PUF-2,  it  is  anticipated  that  the  pressure  prediction  by  PSF 
will  be  better  than  that  by  PUF-2. 
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F  :  Intersection  of  the  profile  base  line  and  the  center  line 
or  generator  line 

P' :  Intersection  of  the  nose-tail  line  and  mid-chord  line, 
or  blade-reference  line 


ITTC  PC 

PUF-2 

Coord. 

oxy 

Pitch 

Angle 

♦p 

Chord 

Length 

c 

gj 

Figure  3  -  Definitions  of  Blade  Surface  Ordinates  Used  in 
ITTC  Propeller  Cooanlttee  and  PUF-2 


SUCTION  SIDE  PRESSURE  SIDE 

Rn=1.1  xIO®  J  =  1.14 


SUCTION  SIDE  PRESSURE  SIDE 

Rn  =  2.6x  106,  J=  1.15 


Figure  6  -  Surface  Flow  Patterns  by  Oil-Film  Test 


TABLE  1  -  PRINCIPAL  CHARACTERISTICS  OF  MODEL  PROPELLER  MP  282 


Diameter  (D)  950.0  mm 

Pitch  (P)  1140.0  mm 

Pitch  Ratio  (P/D)  »t  the  Tip  1.20 

Boss  Ratio  (d/D)  0.198 

Thickness  Ratio  (t/D)  0.05 

Expanded  Area  Ratio  (Ag/AQ)  0.639 

Rake  Angle  0  deg 

Number  of  Blades  (Z)  4 

Blade  Section  MAU 


r/R 

CLE  (Bm) 

CTE  (ani) 

Pitch  Angle  (deg) 

Q 

2 

132.2 

95.8 

62.36 

3 

151.7 

114.6 

51.85 

4 

165.6 

132.8 

43.68 

0. 

5 

174.1 

149.2 

37.38 

0. 

6 

175.3 

164.4 

32.48 

B 

7 

165.6 

175.9 

28.62 

0. 

8 

138.9 

179.6 

25.52 

H 

9 

86.1 

166.2 

23.00 

95 

46.4 

144.2 

21.90 

B 

0 

-59.2 

59.2 

20.91 

p:  blade  width  from  the  centerline  to  the  leading  edge 
p:  blade  width  from  the  centerline  to  the  trailing  edge 
(see  Figure  3) 
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TABLE  2  -  GEOMETRY  OF  MODEL  PROPELLER  MP  282  IN  PUF-2  COORDINATE  SYSTEM 
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TABLE  4  -  TABLE  OF  OFFSETS  OF  MODEL  PROPELLER  MP  282  ABOUT  NOSE-TAIL  LINE 
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TABLE  6  -  COMPARISON  OF  OPEN-WATER  PERFORMANCE  OF 
MODEL  PROPELLER  MP  282 


Experiments 

Numerical 

*T 

10  Kq 

n 

*T 

10  Kq 
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0.9 

- 

- 

- 

0.188 

0.388 

0.702 

0.95 

- 

- 

- 

0.168 

0.348 

0.724 

1.0 

0.1669 

0.340 

0.797 

0.144 

0.312 

0.742 

1.054 

0.1405 

0.294 

0.8025 

0.120 

0.268 

0.753 

1.1 

0.1184 

0.256 

0.803 

0.100 

0.232 

0.754 

1.163 

0.088 

0.205 

0.7865 

0.072 

0.180 

0.729 
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TABLE  8  -  PREDICTED  PRESSURE  DISTRIBUTION  AND  LIFT  COEFFICIENTS  AT  J  =  0.95 
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TABLE  10  -  PREDICTED  PRESSURE  DISTRIBUTION  AND  LIFT  COEFFICIENTS  AT  J  =  1.054 
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TABLE  12  -  PREDICTED  PRESSURE  DISTRIBUTION  AND  LIFT  COEFFICIENTS  AT  J  -  1.163 
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OTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSRDC  REPORTS.  A  FORMAL  SERIES.  CONTAIN  INFORMATION  OF  PERMANENT  TECH¬ 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A  SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY,  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN¬ 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 
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